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ABSTSiCT 


Per  cer-^ain  assumptions,  an  analysis  of  multipla  MAD 
signals  results  in  a  reasonable  estimate  for  the  localiza¬ 
tion  of  a  target  relative  to  the  SAD  platform.  This  is 
achieved  by  using  selective  approximations  to  linearize  an 
initially  nonlinear  problem.  The  simulation  ignores  noise 
and  requires  an  estimate  of  the  lagnitude  of  the  target 
magnetic  moment  components.  Results  indicate  that  the  best 
localization  estimates  are  achieved  when  the  platform  is  on 
cardinal  headings,  and  when  the  target  moment  has  a  strong 
vertical  component. 
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I.  IIIgODOCTION 


A.  BACK6B008D 

Ever  since  the  operational  iatroduct ior.  of  the  MAD 
(Magnetic  Anoealy  Detection)  systea  in  1944,  it  has  been  a 
single  sensor  system  limited  to  indicating  detection  only. 
The  short  ranges  of  the  early  versions  allowed  a  general 
localization  if  a  signal  was  detected.  The  nse  of  multiple 
sensors  should  allow  for  a  good  estimate  of  the  target's 
location  relative  to  the  sensor  platform,  assuming  that  the 
arrangement  of  the  sensors  is  such  that  the  received  signals 
are  linearly  independent  of  one  another.  Previous  attempts 
to  solve  this  problem  have  in  general  focused  on  the  use  of 
only  two  sensors,  but  did  not  develop  a  usable  solution  to 
operational  completion.  »orlc  has  been  dene  on  this  topic  by 
Hynn  et  al. ,  [  Hef .  1],  using  a  superconducting  qradiometer 

array.  The  measured  signals  in  their  work  were  subjected  to 
a  novel  signal  processing  technique  which,  in  the  labora¬ 
tory,  was  used  to  invert  the  dipole  field  equation  to  deter¬ 
mine  both  the  position  and  the  moment  vectors  of  a  dipole 
signal  source.  This  thesis  approaches  the  same  problem  by 
simulating  five  standard  total  field  measurements.  As  shown 
in  Figure  1.1  ,  four  magnetometers  moving  in  parallel  paths 
are  used  to  obtain  these  measurements.  The  advantage  of  this 
approach  compared  to  that  of  Wynn  et  al.  is  that  it  utilizes 
existing  technology  in  conjunction  with  approximations, 
delineated  in  Chapter  3,  in  order  to  produce  a  usable 
result. 
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B.  APPROACH 


The  problem  of  localizing  a  target  from  any  plaxform 
involves  at  least  three  unknowns:  r^  ,  r^  ,  and  r ^  ,  the  compo¬ 
nents  of  the  position  vector  from  the  platform  to  the 
target.  In  the  case  of  using  MAD,  three  additional  unknowns 
m^#  ffly,  and  the  components  of  the  target's  magnetic 

moment,  are  also  involved.  The  result  is  that  there  are  six 
unknowns  that  must  be  solved  for  in  order  to  localirs  a 
target  using  HAD.  Since  a  normal  sized  magnetic  targi  ^an 
be  approximated  by  the  magnetic  field  of  a  dipole  [Hel  ?], 
it  was  decided  to  use  the  dipole  equation  as  dsscrit  in 
[Ref.  3]  to  generate  the  signals  used  in  this  simul- 
Equation  1.1  describes  the  MAD  signal,  and  will  be  deve..jped 
from  the  dipole  equation  in  chapter  2.  Equation  1.1  is  non- 

S{r)*-(m^!/n^y+m^/r3+3(ay:;*-m^+!i^  (r;^^4--yr;^/rs  (eqn  l.l) 

linear  in  the  unknowns.  The  position  vector  in  equation  1.  1 
is  defined  as  going  from  the  target  to  the  platform;  this  is 
a  consequence  of  using  the  dipole  equation,  and  will  be 
explained  later. 

The  cartesian  coordinates  are  defined  in  Figure  1.2  . 
The  z  axis  is  in  the  direction  of  platform  motion,  and  the  y 
axis  is  vertical  down.  The  ccmponents  of  the  earth's  field 
unit  vector,  e,  ,  ,  and  a^  ,  are  assumed  to  be  known  quanti¬ 

ties.  Since  it  is  very  difficult  to  explicitly  solve  a  set 
of  non-linear  equations  in  six  unknowns,  the  problem  will  be 
simplified  by  making  assumptions  and  dividing  it  into  two 
special  cases.  The  first  case  is  a  trivial  on®  in  which  the 
target's  position  is  assumed  known.  For  this  case,  equation 
1.1  is  linear  in  the  three  unknowns  ,  m^ ,  and  m, .  If  the 
three  aeasureients  are  linearly  independent,  then  Cramer's 
Rule,  [Ref.  4],  may  be  used  to  solve  for  these  unknowns.  The 
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second  case  is  the  mere  difficult  one,  and  it  is  th&  prebi-n 
investigated  in  this  thesis.  In  this  case,  the  assumpticr.  is 
made  that  the  target's  magnetic  moment  is  known.  That 
reduces  equation  1.1  into  a  non-linear  equation  in  only 
three  unknowns,  r, ,  r^,  and  r  ^ .  If  rhis  simplified  version 
of  equation  1.1  could  be  manipulated  into  a  linear  form  by 
using  several  sensors  and  appropriate  combinations  of 
signals,  the  position  of  the  target  could  easily  be  solved 
for. 

This  thesis  will  develop  an  approximate  solution  to  this 
problem  by  using  combinations  of  the  total  field  signals  of 
an  idealized  multiple  MAO  system. 


1 1 


Figure  1.1  Spatial  Relationship  of  Sensor  Positions. 
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II.  THBOHI 


The  to-tal  field  magnet oaeter  measures  the  magnitude  of 
the  vector  sum  of  the  aaeient  magnetic  field  oz  the  earth 
and  the  magnetic  field  of  a  magnetic  dipole  moment.  The 
field  of  the  dipole  will  used  to  represent  the  magnetic 
field  of  the  target.  The  magnitude  of  the  vector  sum  is 
called  the  total  field  as  given  by 

K  *  I3e  ♦ 

represents  the  earth's  magnetic  field  vector  and 
represents  the  field  due  to  a  magnetic  dipole.  In  the 
general  case  of  interest  here,  the  magnitude  of  the  dipole 
field  is  much  less  than  the  magnitude  of  the  earth's  field. 
Onier  certain  conditions,  a  useful  approximation  to  the 
total  field  may  be  derived  from  the  geometry  shown  in  Figure 
2.1  . 


Figure  2.1  7ector  Relationship  of  ,  bj^,  and  . 


In  this  figure  a  very  small  vector  is  added  to  a  vary  large 
one.  The  vectors  of  Figure  2,1  are  not  to  scale.  is  to 
be  considered  several  orders  of  magnitude  smaller  than  . 
The  application  of  basic  geometric  formulas  for  the  sides  of 
a  triangle  and  the  angle  between  sides  yields 

"  2bj^BgCOSJif 


for  the  geometry  of  Figure  2.1  .  la  terms  cf  0  this  equa-.ioa 
becomes 

3^2  =  3,2  ♦  ♦  2bj^B^cose 

The  3^2  term  can  be  factored  out  of  the  right  hand  side  of 
the  equation  to  produce 

♦  (bA/Bg)2  «■  2(b^/Be)cos9} 

Talcing  the  square  root  yields  equation  2.1  . 

(bji/Be)  *  +  2(b//Be)cos9  3Vfe)  (egn  2.1) 

Since  by  assumption  b^  is  much  less  than  3, ,  the  oerm 
(b^/B,)  2  is  much  less  than  2  (bj/Be)  cos9  and  can  be  cor.sii- 
ered  negligible.  This  reduces  equation  2.1  to 

at  =  B,  {CU2(bj^/Be)  C0S9]V<!}  Oqn  2.2) 

The  general  expansion  formula  given  by 

(l+x)*'  =  1  ♦  nx  ♦  n(n-1)x2/2!  +  Higher  Order  Terms 

can  be  applied  to  equation  2.1  .  For  the  case  here,  a=1/2, 
and  x=  (bj^/Bg  )2 >2  (bj_/B^)  C039  ,  Substituting  these  values  into 
the  general  expansion  formula  gives  equation  2. 1  as 

3^  =  B,  {1  ♦  .5C(bj^/B,)2  ♦  2(bi/Be)  COS9] 

-  1/8[  (bj^/B*  )  2  ♦  (2bji/Be)  C0  39  ]2  ♦  H.O.T.) 

Squaring  the  appropriate  term  and  multiplying  through  by  the 
coefficients  results  in 

3^  =  B,  {1  ♦  .5C  (bj^/Bj)  2]  ♦  (b^/BelC039  -  1/3[b^/Be]^ 

-  4/8(  bj|_/B^  ]3C0S9  -  4/8[b4_/B,  ]2cos29  ♦  H.  O.T.} 

The  cubed  and  higher  order  terms  will  be  considered  to  be 
negligible.  Dropping  those  terms  and  rearranging  the 
remaining  ones  yields  the  form 
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=  Be  {1  ♦  (bj^/Be)cose  -  . 5  (b^/Be )  ^cos^e  ♦  .5(b^/5e)2} 

Using  the  relation  3in*9  ♦  cos*©  =  1,  this  radaces  tc 

®t  “  ®e  ♦  (bj[/Be)cos9  *■  . 5  (bj^/Be )  *sin*©) 

After  multiplying  through  by  3* »  the  second  term  car.  be 
written  in  tens  of  the  dot  product  bj^»Be/3e  .  This  is 
reduced  to  rf  the  unit  vector  "t  is  defined  as 

and  the  resulting  equation  is 

3^  =  B,  ♦  b^«t  *■  .  5  (b^/Be )  2sin*© 

Since  bj;_  is  on  the  order  of  InT  while  3^  is  on  the  order  of 
50,000  nT,  the  error  term  involving  the  sine  is  negligible. 
For  the  stated  assumptions,  equation  2.1  can  be  repres#»nted 
by 

(eqn  2.3) 

In  this  thesis,  3|,  will  be  assumed  constant.  In  addition, 
this  simulation  will  assume  that  no  noise  exists.  Thus,  the 
total  field  magnetometer  will  measure  the  anomaly  in  the 
earth's  field  defined  by 

S  (r)  =5 

In  this  expression,  bj^,  represents  the  magnetic  field  of  a 
dipole.  is  defined  £Bef.  3]  by  the  equation 

b^  <?)  =  -mA^  ♦  3(m«r)r/r5  (aqn  2.4) 

m  is  the  magnetic  dipole  moment  in  units  of  nTft^.  r  is  the 
position  vector  to  the  field  point  and  is  measured  in  feet. 
If  b^  is  much  less  than  B«  ,  than  aquation  2.3  and  equation 
2.4  result  in  the  aquation 
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S  (r)  =-  (n»e) /r  ^  +3  (m •r)  (f  •a )  /r * 


(eqn  2.5) 


This  equation  is  an  alternatire  basis  for  the  derivation  of 
the  Anderson's  Functions,  found  in  [Bef.  5].  The  Anderson 
Functions  are  a  coamonly  used  approxiaation  for  the  signal 
received  by  a  SAD  sensor.  The  coiputer  simulation  in  this 
thesis  is  based  on  equation  2.5  directly,  since  it  is  in  a 
representation  which  is  easier  to  use  than  the  Anderson's 
function  representation  of  a  SAD  signal. 
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III.  HATOSJ  Of  THE  PHOfiiEH 


k,-. 


This  thesis  Kill  develop  and  demonstrate  an  approxima¬ 
tion  that  generates  a  reasonable  solation  of  the  target's 
position  relative  to  a  platform  bearing  several  "sensors”  to 
"measare"  the  nagnetic  signal  as  the  platform  makes  a 
straight  line  encoanter  near  the  target.  The  assumptions  are 
as  follows: 

1. )  The  dipole  equation  is  a  valid  representation  of 
the  magnetic  field  of  the  target. 

2. )  bjL  is  several  orders  of  magnitude  smaller  than  3,. 

3. )  Heasur aments  are  done  using  only  straight  line  en¬ 
counters. 

4. )  The  target's  magnetic  dipole  moment  is  a  known 
vector. 

5. )  The  earth  magnetic  field  vector  is  known. 

6. )  The  gradient  of  the  earth's  field  is  zero. 

7. )  No  magnetic  noise  of  any  kind  (environment,  sensor, 
etc)  3xists. 

8. )  The  measurements  are  oriented  along  four  lines 
which  are  oriented  with  respect  to  the  flight  path  as 
shown  in  Figure  1.1  . 

Since  the  dipole  equation  is  written  with  the  moment  at  the 
origin  of  a  coordinate  system,  the  calculations  of  the 
signal  field  values  used  in  the  simula-ion  are  made  wich  the 
target  located  at  the  origin.  &s  will  be  ?een,  the  program 
performs  a  coordinate  transformation  in  order  to  produce  the 
output  of  target  position  relative  to  the  platform.  After 
this  coordinate  transformation,  the  origin  is  located  at  the 
sensor  platform  and  moves  with  it. 

The  five  measurements  are  made  along  four  separate 
flight  paths,  so  only  four  different  "sensors”  are  involved. 


18 


Sensor  4  could  be  replaced  by  inducing  a  time  delay  in  the 
measurement  of  sensor  1.  The  signals  generated  at  these 
measurement  points  are  measured  at  each  time  step.  The  posi¬ 
tion  vectors  can  all  be  described  relative  to  the  origin, 
located  at  the  target,  in  terms  of  the  position  vector  of 
sensor  1.  The  measured  signals  are  represented  by  the 
follouing  equations: 

S  (f1 )  =- {  (mj^+a^+ n^)  (r1)  2-3  ( ( xa|,+  yy  } /r  1  s 
S(r5)=-{(m^+a,ym^  (r2)  2-3  (y^i^y-d]*^^)  (xy[  y-d  ]yz^  ) /r  2« 

S  (r5)=- { (nyyaye^mi^  (r3)  2-3  ( mj[x-d  ]+y +  11^)  ([  x-d  ]e|,+  yy  z^  } /r3s 
S(r1l)=-{(a^ii^yn^  (r4)  2-3  (nyt+iy+i^  z-d  ])  (x9/yy[  z-d  } /r  4* 
S(r^)=-{(ayayymj^  (rS)  2-3  (;^x+d  +  ([  x+d  le  +  yey+z^  ) /r 55 

For  convenience,  the  following  dummy  variables  are  used  in 
the  computer  program  and  are  defined  as: 

G  =  a*s^  +  myey  ♦m^a^ 

71  =  3m^e,t  -  G 

72  «  3  (m,«y  ♦  ) 

73  a  3my€y  -  G 
74=3  (m^ej  ♦  m^Sy  ) 

75  =  3m,e^  -  G 

76  =  3(m,e^  -  m^e^ ) 

51  ^  S(f^) 

52  =  S (f^) 

53  a  S(?3) 

54  a  S  (fft) 

35  a  S(r§) 

The  use  of  these  variables  reduces  the  above  equations  to 
the  following  form,  which  are  the  defining  equations  for  the 
signals  in  the  simulation. 

S1=  (71x2+7  2xy+7  3y2  +  74yz+75  z2+76xzj  /r  1  = 
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S2=  (V1x2+V2x7-72xa  +  V3y*-2V3yd+73d2  +  V47z-V4dz+V5z2  +  V6x2) /r2s 
S3=  (V1x2-2V1xd+71d2  4-V2xy-?2dy+73y2«-y4yz+V522+76xz-76az)/r3  5 
S4=  (V1x2+7  2xy+73y2  +  74yz-74yd+75z2-2752a-*-75d2  +75x2-76x3)  /r45 
S5=  (71x2 +271x3  +  71  d2+72xy+7  2dy +7372 +  7 4y z+7522 +76 xz+V6.iz)  /rS* 


The  unknowns  of  interest  are  x,  y,  and  z,  which  are  related 
to  the  position  vectors  as  follows: 


rl 

a  Xt  + 

yt 

+  z1c 

r1 

II 

♦ 

(y- 

d)T  ♦ 

zlc 

r5 

a  (X-d)t  + 

yt  ♦ 

zlc 

n 

a  X”?  + 

y^ 

+  (z- 

d)1c 

r$ 

=  (x+d) 

t  + 

yf  * 

zlc 

(eqn  3.1) 
(eqn  3.2) 
(agn  3.3) 
(eqn  3.4) 
(eqn  3.5) 


To  linearize  the  problea  it  will  be  assumed  tha*  the 
denominators  of  the  equations  for  Si,  S2,  S3,  34,  and  S5 

given  above  are  all  equal  to  r,  i.e. , 
r  =  r1  =  r2  *  r3  =  r4  =  r5 

Tha  equations  for  Si,  S2,  S3,  34,  and  S5  with  the  denomina¬ 
tors  equal  can  then  be  subtracted  from  one  another  to  yield 
the  following: 


S2-S3={(72-271)  x+  (273-72)y+  (74-76)  z+  (71-73)d)  (-d/r®) 

52- S5={(V2+271)  x+  (273+V2)y  +  (74+75)  z+  (71-73)  dj  (-d/r*) 

53- S1=  {271x  +  72y+76z-71d)  (-d/r*) 

S5-S1=  (- 271  X- 72 y- 762-713)  (-d/rS) 

S3-S4={  (271-7  6)  x+  (72-74)  y+ (76-275)  z+  (75-71)3)  (-d/rS) 
S5-S4=  (-  (271  +  76)  x-(V2+74)  y-  (76  +  275)  z+  (75-71)  d)  (-d/r*) 

Since  tha  denominators  are  assumed  to  ba  identical,  ii  these 
six  equations  are  divided  in  pairs  the  coefficients 
involving  the  non-linear  r*  term  and  the  quantity  d  are 
canceled  out.  This  results  in  tha  following  equations  that 
define  the  new  variables  A1  ,  A2,  and  A3  as 


A1  *  (S2-S3)  /  (S2-S5) 
A2  a  (S3-S1)  /  (S5-S  1) 
A3  a  (S3-S4)  /(S5-S4) 


(eqn  3.5) 
(eqn  3.7) 
(eqn  3.3) 


Daficing  the  fallawing  variablas  for  coa^eniance. 

Cl  =  A1(V2+2V1)  -  (V2-2V1) 

C2  =  A1(2V3+72)  -  (2V3-72) 

C3  =  A1  (74+76)  -  (74-76) 

K1  =  [  (71-73)  -  A1  (71-73)  ]d 

C4  =  A2  (271)  ♦  271 

C5  =  A 2  (72)  ♦  72 

C6  =  A2(76)  +  76 

K2  *  [  71  -  A2  (71)  ]d 

C7  s  A3  (76+271)  -  (76-271) 

C8  =  A3  (74+72)  -  (74-72) 

C9  =  A3  (275+76)  -  (275-76) 

K3  =  [  (71-75)  -  A3  (71-75)  ]d 

equaticas  3.6  -3.9  can  be  writ-en  in  zhe  following  form  as 

Clx  +  C2y  +  C3z  =  K1  (eqn  3.9) 

C4x  +  C5y  +  C6z  =  K2  (eqn  3.  10) 

Clx  *  C8y  +  C9z  «  K3  (eqn  3.11) 

These  three  equations  ara  linear  in  terms  of  x,  y,  and  z.  In 
addition,  they  ara  linearly  indepanient  as  shown  by  using 
the  program  in  Appendix  B.  A  sample  ounpun  demonstrating  the 
linear  independence  is  shown  in  Appendix  C.  Since  they  are 
linearly  independent,  Cram3r*5  Rula  can  be  used  to  solve 
equations  3.9  -  3.11  for  x,  y,  and  z.  The  computer  program 
in  Appendix  A  was  written  in  terms  of  ohese  equations.  The 
program  assumes  that  the  magnitude  of  the  target's  magnetic 
moment  is  5x10*  nlft^.  A  r eprasentat iva  moment  was  selected 
by  using  the  values  listad  by  Fromm  in  [Ref.  2],  who  lists  a 
moment  of  10*  to  2x10*  cgs  units  far  a  submarine.  One  cqs 
unit  is  approximately  3.35  nlft^.  The  value  used  in  the 
program  is  a  rounded  out  average  of  Fromm's  values. 
Computations  are  done  using  the  standard  convention  with  the 
target  dipole  at  the  origin.  Howevar,  as  indicated  earlier. 
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;  tha  output  is  corrected  by  a  coordinate  transf oranati cn  sc 

'C^  that  the  platfora  is  at  the  origin,  and  the  target  position 

is  given  with  respect  to  the  platform. 

k 


I?.  DISCaSSIOH 


A.  GEHSBAL 

The  program  in  Appendix  B  was  written  to  check  whether 
the  results  generated  by  the  approximation  equations  were 
reasonably  close  to  the  actual  coordinates  of  the  target. 
The  program  was  modified  into  the  fora  shewn  in  Appendix  A 
to  graphically  provide  a  comparison  between  the  calculated 
and  the  actual  position  coordinates.  Initial  parameters  for 
the  knowns  in  the  equations  were  a  platform  heading  of  30° 
magnetic,  an  earth  field  vector  of  TO®  down  from  the  hori¬ 
zontal,  and  a  target  dipole  moment  vector  with  a  vertical 
component  50®  down  from  the  horizontal  and  a  horizontal 
component  oriented  355®  from  magnetic  north.  These  initial 
parameters  have  no  special  significance,  and  merely  repre¬ 
sent  a  convenient  starting  point.  As  discussed  earlier,  the 
magnitude  of  the  target  dipole  moment  was  set  at  the  value 
of  5X10®  nTft®.  The  spacing  between  measurement  positions, 
d,  was  taken  as  50  feet,  and  the  platform  passed  directly 
overhead  the  target  at  1000  feet.  This  and  all  subsequent 
runs  were  started  at  5000  feet  prior  to  CPA  (closest  point 
of  approach)  with  the  platform  advancing  50  feet  per  time 
step.  One  position  calculation  was  completed  at  each  time 
step,  the  simulation  ceasing  when  the  platform  was  5000  feet 
past  CPA.  For  interpretation  purposes,  calculated  target 
positions  are  shown  as  the  sensor  platform  moves  10,000  feet 
past  the  target  in  a  straight  line  encounter,  with  CPA 
occurring  at  the  5000  foot  point.  Neither  the  CPA  nor  the 
targ9t*s  location  are  in  any  way  considered  known  in  the 
localization  process.  However,  as  the  area  of  interest  is 
the  behavior  of  the  simulation  close  to  CPA,  the  closest 
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Figure  4,  1  Initial  Bun,  x  position. 

nt  of  approach  is  always  centersd  in  tha  follcwina 
ures. 
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Figure  4.3  Initial  Hun,  z  position. 

this  thesis  investigates  straight  line  encounters  only, 
x(actual)  and  y  (actual)  cocrdinates  will  remain  cons- 
throuahout  any  run  for  the  coordinate  oriantaticn  used. 


x(actual)  would  be  constant  at  zaro,  y  (actual)  would  be 
constant  at  1300  feet,  while  z  (actual)  would  progress  frcm 
5000  feat  to  -5000  feat.  Figures  4.1-  4.3  show  a  compar¬ 
ison  of  the  actual  and  calculated  target  positions  as  the 
platfora  completes  tha  initial  run  described  earlier.  It 
will  be  noted  that  the  calculated  and  actual  x,  y,  and  z 
positions  are  almost  identical  at  the  start  of  the  run.  This 
is  a  coincidence.  "Crossover”  points  such  as  this  occur  for 
all  three  coordinates  simultaneously  during  a  run,  falling 
closer  or  further  away  from  CPA,  depending  upon  the  choice 
of  parameters. 

In  this  ininial  run,  it  is  obvious  that  the  "worst" 
estimates  occur  for  the  y  coordinate.  This  is  less  signifi¬ 
cant  than  it  seems  at  first  glance,  however,  as  the  y  coor¬ 
dinate  is  the  vertical  separation  between  the  platform  and 
the  target.  If  the  target  is  submerged,  this  value  is 
unknown,  but  the  altitude  of  the  platform  is  known.  The 
target  can  be  no  closer  to  the  platform  than  this  altitude, 
since  the  target  is  physically  unable  to  rise  above  the 
ocean  surface.  As  an  example,  if  a  target  depth  of  200  feet 
is  assumed,  the  initial  y  separation  of  1000  feet  forces  tha 
altitude  in  this  simulation  to  be  830  feet.  Osing  this  value 
as  a  bound  on  y ,  it  can  be  seen  in  Figure  4.4  that  the  y 
calculated  position  is  somewhat  improved.  Accordingly,  it 
seems  that  the  particular  arrangements  of  sensors  shown  in 
Figure  1.1  does  show  promise  in  the  problem  of  utilizing 
this  geometry  of  Multiple  MAD  for  accurate  localization 
purposes.  Accordingly,  variations  of  parameters  were  tried 
in  order  to  determine  what  effect  changing  them  would  have. 
In  all  cases,  parameters  not  specifically  mentioned  in  a 
section  will  retain  the  values  as  for  the  initial  run 
described  above.  The  singularity-like  behavior  occuring  at 
various  places  on  the  graphs  are  a  result  of  the  determinant 
of  the  denominator  in  the  Cramer's  Rule  subroutine  having  a 
value  close  to  zero. 
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B.  HEADING  CHANGES 


In  ”his  saction,  the  localization  process  will  be  ir.ves- 
tigated  as  a  function  of  heading.  Compared  to  the  init-al 
conditions  described  above,  the  crossover  points  were  found 
to  be  moved  further  out  past  5003  feat  prior  to  CPA  on  a 
heading  of  Oo.  While  a  specific  example  of  this  is  not 
shown,  the  behavior  of  tha  crossover  points  for  the  platform 
headings  of  33o  and  60o  may  be  observed  by  comparing  Figures 
4. 1  -  4.4  to  the  respective  portions  of  Figura  4,5  .  The  x 
position  estimates  were  better  for  a  North -South  heading 
than  for  any  other  platform  heading,  although  rhe  estimates 
on  East-West  headings  wera  also  very  good.  For  a  90°  heading 
thare  are  two  crossover  points  locatad  equidistant  at  1200 
feet  from  CPA.  Calculated  positions  on  all  headings  are 
generally  good,  excepting  those  around  6Co,  where  the  x 
estimates  shows  a  large  discrepancy  at  distances  greater 
than  3000  feet  prior  to  CPA.  The  y  astimates  are  not  signif¬ 
icantly  affectad,  but  the  z  estimates  are  good  only  within  ± 
1200  faet  of  CPA.  Nevertheless,  avan  in  this  worst-case 
situation,  reasonable  calculated  coordinates  are  available 
for  localization  along  a  significant  portion  of  the  track. 

As  shown  in  Figures  4.5  -  4.8,  headings  in  the  quadrant 
go  -  90  0  reasonably  describe  what  happens  on  any  heading,  as 
the  graphs  for  analogous  headings  (53o,  120o,  240°,  and  300o 
in  this  example)  have  similar  shapes,  although  the  crier.ta- 
ticn  may  differ  between  quadrants.  These  similarities  are 
apparently  due  to  the  fact  that  analogous  heading  in  quad¬ 
rants  1  and  3  and  quadrants  2  and  4  are  in  fact  reciprocal 
headings,  and  so  one  would  expect  the  graphs  to  maintain  the 
same  shape;  however,  the  ending  point  of  one  would  represent 
the  starting  point  of  the  other.  The  relationship  between 
adjacent  quadrants  probably  exists  for  similar  reasons. 
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C.  TARGET  DIPOLE  MOMENT  CHANGES 


In  this  section,  the  horizontal  component  of  m  vas 
oriented  in  the  Qo  (magnetic  North)  direction  vice  3550,  and 
the  vertical  component  of  m  was  varied.  The  best  positioning 
was  obtained  for  m  orien^-ed  vertically  up,i.e.,  no  hori¬ 
zontal  component,  and  an  example  for  that  case  is  shown  in 
figure  4.9  .  Poor  positioning  was  obtained  for  the  vertical 
component  of  ■  pointing  between  the  horizontal  and  45®  up. 
The  worst  case  shown  is  Figure  4.10,  whers  z  is  oriented 
horizontally  pointing  due  north.  Bith  the  moment  oriented  at 
an  angla  below  the  horizontal,  reasonable  positioning  esti¬ 
mates  ware  obtained.  There  was  less  variability  in  posi¬ 
tioning  for  moments  pointing  below  the  horizontal  than  for 
those  pointing  above  the  horizontal. 

The  results  for  variations  in  the  orientation  of  the 
target  dipole  moment  suggests  that  the  positioning  estimates 
are  less  sensitive  to  moment  orientation  changes  in  the 
horizontal  plane  as  compared  to  the  vertical  plane.  For 
example,  by  keeping  the  vertical  component  of  m  constant  at 
an  angle  cf  50®  down  while  varying  the  horizontal  component 
direction  it  was  fcund  that  the  horizontal  component 
affected  the  positioning  much  less  than  variations  in  the 
vertical.  Illustrated  in  Figure  4.11,  the  x  estimates  of 
position  improve  as  the  horizontal  component  approaches  the 
platform  heading,  while  the  y  estimates  improve  as  the  hori¬ 
zontal  component  approaches  alignment  with  magnetic  North. 
The  latter  case  is  not  shown. 


D.  CHANGES  IN  SENSOB  SPACING 

Changing  the  inter-sensor  spacing  had  at  most  a  wealc 
effect  upon  the  position  output,  even  with  spacing  as  short 
as  5  feet  or  as  long  as  200  feet.  This  would  be  expected  in 
a  computer  simulation  where  sufficient  precision  is 
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available.  For  a  spacing  as  gr=at  as  1030  fast  sdh?  afirC* 
bacama  noticabl*.  in  nh®  raal  i/orld,  =ansor  sparine  W21il 
inflaanca  '■h?  localization  “stiaiatss. 

E.  ALTITODE  CHANGES 

For  a  straight  cverhaad  pass,  positioning  information 
appears  best  if  a  pass  with  a  500  to  1000  foot  separation 
between  the  targat  and  the  platform  is  made  at  CPA.  At 
greater  altitudes,  the  accuracy  suffers  slightly  with  the 
increase  in  altitude,  while  at  lessar  ranges,  usable  posi¬ 
tioning  information  is  obtained  only  relatively  close  to 
CPA.  For  example,  for  a  vertical  separation  of  400  feet, 
positioning  data  was  very  inaccurate  at  slant  ranges  greater 
than  1000  feet. 

F.  NOH-OfEBHB&O  PASSES 

This  program  gives  poor  results-  for  passes  that  are  nc^ 
directly  over  or  directly  off  to  one  side  of  the  target.  The 
worst  situation  observed  in  any  variation  is  illustrated  in 
Figure  4.12  .  To  generate  this  figure,  a  pass  was  made  such 
that  the  actual  x  coordinate  was  approximately  equal  ■^.c  the 
actual  y  coordinate.  The  best  positioning  was  achieved  by 
overhead  passes,  but  low  off  to  the  side  passes  can  yield 
good  results  as  well. 

G.  AVERAGES 

The  averages  resulting  from  a  number  of  runs  are 
presented  in  Table  I  .  These  selected  averages,  covering 
all  of  the  variations  discussed  in  this  chapter,  show  little 
consistency.  This  is  probably  due  in  part  to  the  existence 
of  the  singularities.  The  variability  is  also  a  function  of 
the  limits  of  travel  for  the  runs.  Nevertheless,  for  most 
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cases,  they  d3  shdw  reasonable  values,  and  by  themselves  io 
provide  limited  localization  data.  The  problem  for  inflight 
analysis  would  be  to  determine  the  bounds  of  a  run  in  zhs 
absence  of  any  other  target  positioning  data. 

B.  IHPOBTAHCB  OF  CORSISTEHCT 

in  important  observation  resulting  from  the  analysis  of 
the  graphs,  including  ^hose  not  shown  as  examples  in  this 
chapter,  is  than  the  rate  of  change  of  the  calculated  coor¬ 
dinates  gives  some  indication  of  the  validity  of  the  posi¬ 
tioning  data.  If  the  x  and  y  graphs  have  a  slope  close  to 
zero,  and  if  the  z  graph  has  a  slope  of  -45®  for  the  calcu¬ 
lated  coordinate,  then  positioning  data  is  likely  to  be 
comparatively  accurate.  If  all  three  graphs  fulfill  this 
condition  simultaneously,  then  positioning  data  appears  to 
be  usable. 

1.  IflPLBBEHTATION 

The  possible  implementation  of  the  process  described  in 
this  thesis  would  be  to  assume  a  magnetic  moment  for  the 
target,  then  run  the  signal  data  through  the  processing  to 
generate  positioning  data,  and  then  compare  these  results  to 
those  obtained  by  using  the  assumed  moment  to  generate  a 
signal  via  this  simulation.  If  the  agreement  between  results 
were  good,  cne  could  assume  that  one  had  assumed  the  correct 
target  moment;  if  not,  one  would  assume  a  different  moment, 
and  go  through  the  process  again,  until  one  had  a  reasonable 
target  moment,  which  would  be  valid  at  least  until  the 
target  changed  course.  Osing  a  computer  these  comparisons 
could  be  made  reasonably  quickly,  so  this  process  might  have 
operational  possibilities. 


4  1 


?.  COICLOSIOB 


For  the  assumpticns  mads,  the  positioning  ^stimatas  are 
reasonable.  &lso,  the  a/erages  talcen  of  the  three  coordi¬ 
nates  as  calculated  throughout  the  various  runs  show  rhat 
the  X  calculation  averages  are  in  general  ar  worst  wirhin 
500  feet  of  the  actual  x  value,  while  the  worst  y  calcula¬ 
tion  averages  tend  to  be  within  10  09  feet  of  the  actual  y 
value.  The  bast  calculated  x  and  y  values  were  within  100 
and  200  feet,  respectively,  of  the  true  values. 

Trends  observed  were  as  follows: 

1.  It  is  best  to  fly  the  sensor  platform  directly  over 
the  target,  or  alternatively  to  have  a  'distant* 
(between  1000  and  2000  foot)  CPA  off  to  one  side  while 
minimizing  the  target  to  platform  separation  in  the  y 
direction.  The  worst  situation  is  when  the  target  is  at 
a  point  that,  at  CPA,  has  a  depression  angle  of  45o  with 
respect  to  the  platform. 

2.  Best  vertical  separations  are  those  that  will  yield  a 
500  tc  1000  foot  vertical  distance  between  target  an  i 
platform.  Short  vertical  separations  should  be  avoided 
if  good  localization  is  desired. 

3.  Best  results  are  obtained  if  headings  are  northerly 

or  southerly.  Poor  results  are  achieved  wi^h  headings 
close  to  600,  1200  ,  24 Oo,  and  3030.  East-West  headings, 

while  not  as  good  as  North-South,  do  provide  good  posi¬ 
tioning  data  close  to  CPA. 

4.  Good  results  may  be  expected  if  T.4.>.r9  is  a  large 
vertical  moment,  especially  if  it  is  vertically  up. 
Poor  positioning  results  are  obtained  with  this  simula¬ 
tion  if  the  moment  has  a  large  horizontal  component. 


5.  Results  are  slightly  better  if  the  target  aomert  is 
lined  up  with  either  the  axis  ef  flight  or  wi*h  the 
earth's  field,  although  this  effect  is  often  dominated 
by  the  more  gross  effects  caused  by  the  ether  variables. 

A  suggested  approach  for  further  study  is  to  inject 
noise  into  the  simulation  in  order  to  see  how  robust  it  is 
in  the  accuracy  of  the  position  calculations.  If  the 
results  remain  reasonable,  this  approximation  or  an  altered 
form  of  it  may  one  day  have  real-world  applicability. 
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APPENDIX  A 


COMPUTER  PROGRAM  FOR  PLOT  OUTPUT 
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COMPUTER  PROGRAM  FOR  NUMERICAL  OUTPUT 
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4E2M2IX  c 

BX&HSLB  OOTPOT  0?  &PPEBDIX  B 


Aa  axaipla  of  the  output  of  taa  program  ia  Appar.dix  B  is 
shown  below.  The  toral  output  for  a  single  time  step  is 
prasented. 


SENSOR 

NUMBER 

1 

2 

3 

4 

5 


EQNS 
USED 
12  3 


SLANT 

RANGE 

5G50.0 

506C.1 

5050.2 

5099.0 

505Q.2 


ACTUAL 


ScNScO 
SIGNAL 
0.60540“ 04 
0.1  5960“  C3 
0.2682  0“C4 
0.3940D“C4 
0.93660-C4 


CALCULAT  EO 
SIGNAL 
0.6054E“04 
0.1612£“03 
0.2683E“04 
0.4135  E“04 
0.9  36dE“04 


CL AT  ION 

TRUE 

APPROXIMATE 

S2-S3 

0.1327C-03 

0.i345E“03 

S2-S5 

0.6590C-04 

0.6748E-04 

S3- SI 

-.3372C-04 

“.3371E-04 

S5-S1 

0.3312C-04 

0.331 4E-04 

S3-S4 

“.1258C-04 

“.i452E“04 

S5-S4 

0.5426C-04 

0.5233E-04 

X 

POS 

0.0 


y 

POS 

>1000.00 


Z 

POS 

-4950.  00 
jr42Ai*.Si 


The  first  fi/e  rows  of  data  above  are  for  the  appropriate 
••sensor”  as  indicated,  the  slant  range  from  that  sensor  to 
the  target,  the  received  signal  as  generated  by  equation  2.5 
and  the  signal  that  would  be  generated  using  the  approxima¬ 
tion  described  in  Chapter  3  to  linearize  the  problem.  The 
next  six  rows  show  the  results  of  taking  the  differences 
between  the  received  and  calculated  (for  comparison  purposes 
only)  signals,  listed  under  the  headings  of  TF.OB  and 
APPSOXISATS,  respectively.  The  final  two  lines  show  the 
actual  coordinates  of  the  platform  relative  to  the  target 
(the  coordinate  transformation  has  not  been  done  at  this 
point)  as  ccmpared  to  the  calculated  position.  These  are  the 
points  used  by  the  program  in  Appendix  A  to  plot  the  graphs 
after  the  coordinate  transformation  is  completed. 
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Th®  second  section  of  the  output  is  shown  below 


JATI  jli 

1 

2 
3 


CJN3TANT 

U*5  J0«' J3 
'~Q«33u+03 
-J*&o0*02 


X 

COtFFiCI£NT 
-0. 300+J2 

Oa 

-Oa  iiD-*-0i 


y 

CQEFFICIdNT 
“Oa  9cU*0  1 

Ja  l^U*-QQ 

“Oa  + 


-  i.  . 

C  Cc  F  F  X  C  a  c  <\  T 
0  a J i 

Wa 

Ua  33u4- 01 


^ATIQ 


1/2 

“Oa  1  3  fc+0  1 

“Oa  22E+03 

“Oa  7Lc* 02 

Oa  Ooc**  02 

1/3 

“0a58c+01 

Oa lAc^Oi 

Oa‘t'siC  +  0  J 

0  a  pOC-^OO 

2/3 

Oa  44c+  01 

“Oa  o3c“0a 

“Oa  6'3»C“02 

Oa 13c“U 1 

The  first  three  lines  of  the  output  represent  the  coeffi¬ 
cients  of  the  equations  3.6  -  3.8  in  the  form  of 

CONSTANT  =  X  COEPF(X)  +  Y  COEFF(y)  *■  Z  COEFF(z).  In  order  to 
determine  the  linear  independence  of  these  equations,  the 
ratios  of  the  coefficients  are  taXen  in  pairs.  Since  it  may 
be  seen  from  the  last  three  lines  of  the  output  that  the 
ratios  of  the  constants  and  coefficients  differ  signifi¬ 
cantly  fcr  all  three  equations,  it  ran  be  concluded  that  the 
equations  are  linearly  independent.  This  holds  true 
throughout  any  given  run,  which  may  be  verified  by  printing 
out  all  time  steps  for  that  run. 
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iPPBHDlX  D 
CSAilEB'S  ROLE 

A.  GEREBAL 

Cramer's  Rule,  a  mathod  for  solving  systems  of  linear 
equations  in  t»o  or  more  unknowns,  states  that  each  unknown 
can  be  expressed  as  the  ratio  of  two  determinates.  The 
general  form  of  a  solution  is  shown  below. 

Ax+By+Cz=0 
Ex+Py+G2=H 
I  x+Jy+Kz=L 

I  A  B  C 

D  =  I  £  F  G 

I  J  K 

0  B  C 

H  F  G  I  /  D 

L  J  K  I 

A  0  C 

y  =  I  E  H  G  /  D 
I  I  L  K 

ABO 

Z  =  £  F  H  /  D 

I  J  L 

B.  EZAEPLE 

The  following  example  was  taken  from  [Ref.  4].  Solve; 
3x-2y+2z*7 
x+  y+  z*6 
2x-  y-2z*2 
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Let  D  be  the  5enoiinator 


D 


X 


7 


z 


13-2  2  1 

I  1  1  I 

I  2-1-2  1 

17-2  2  1 

6  11 
2-1-2 


3  7  2 
1  6  1 
2  2-2 


3-2  7 

1  1  6 
2-12 


-17 

-51 

-17 

-3U 

-17 

-17 

-17 


3 


2 


1 


CHECK: 


3x-2y+2z*7 
3(3)-2<2)t2(1)  =7 
9-at2=7 
7  =  7 


xty+2=6 
3+2+1+6 
6  =  6 


2x-y-2z=2 

2(3)-2-2(1)=2 

6-2-2=2 

2=2 


C.  SOBBOUTINB  "CBAnEfi*S  BOLE"  CHECK 

Runaicg  the  values  in  Section  3  through  the  subroutines 
of  the  programs  in  Appendix  A  and  Appendix  B  yielded  the 
same  results  as  the  example. 
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